A way to solve BCS-type pairing model 
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We propose a way to solve BCS-type pairing model by to exactly solve its spin-analogy in the sub- 
space. The advantages of our method are to avoid to directly deal with the approximate procedure 
and to transfer an exponentially complicated problem to a polynomial problem. Moreover, it builds 
the theoretical foundation of quantum simulation of pairing model and can be used to check the 
precision of quantum simulation. It is also helpful to understand the many-body quantum theory. 
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Recently, the Bardeen, Cooper and Schrieffer (BCS) model for superconductivity .1] has been connected with the 
problems in different areas of physics such as superconductivity, nuclear physics, physics of ultrasmall metallic grains 
and color superconductivity in quantum chromodynamics. It is well-known that BCS-type Hamiltonian is able to be 
studied by its spin-analogy form 00. I n fact, a famous difficulty of exact solution of the spin-analogy of pairing 
model is, in the large-N limit, the obsession of the exponentially complicated problem. In addition, the relation the 
exact solution of spin-analogy of BCS-type Hamiltonian with the spectrum of BCS-type Hamiltonian still have to be 
correctly understood. In this letter, we propose a way to solve BCS-type pairing model by to exactly solve its spin- 
analogy in the subspace. After building the clear relation between the approximate form of BCS-type pairing model 
and the subspace of the spin-analogy of pairing model and finding out the general form of BCS-type Hamiltonian in 
this subspace, we are able to obtain its exact solution in the subspace, that is the eigenvalues of BCS-type pairing 
model under the known physical approximation. The advantages of our method are to avoid to directly deal with 
the approximate procedure and to transfer an exponentially complicated problem to a polynomial problem. It is also 
helpful to understand the many-body quantum theory. 

Just well-known, a possible way to overcome the exponentially complicated problem is to use the quantum simulation 
0. L.-A. Wu et al. reported a NMR experiment scheme performing a polynomial-time simulation of pairing 
models. We also give out a proposal of quantum simulation 7| which can obtain the differences of energy levels of 
pairing model. Moreover, we have implemented our proposal for the simplest system of two qubits in experiment. 
It must be emphasized that one has to back (map) to the simulated physical system 0, that is, one would like to 
know how the spectrum obtained from quantum simulation corresponds to the physical spectrum in the simulated 
physical system, and then one can arrive at the final aim of quantum simulation. Based on such a cause, this letter is 
\ also be used to shows what is the theoretical foundation of the quantum simulations including our pro pos al and can 

be used to check the precision of quantum simulation as well as the physical significance of the results pH | . 
' Let us consider the reduced BCS Hamiltonian 0,13: 
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HBCS = "' O F ( Um + ~ V Y C ™C-m c -lCl (1) 
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where n± m = c± m c± m are the electron number operat ors, ci (c m ) is the fermionic creation (annihilation) operator. 
The coupling coefficient is simplified as a constant V [ill Il2j . Note that the summation indexes m = 1, 2, ■ ■ • , N 
represent all of relevant quantum numbers, and the electron pairs are labelled by the the quantum number m and 
O ■ —to, according to the Cooper pair situation where the paired electrons have equal energies but opposite momenta and 
spins: m = ( k , |) and —to = (— k , j). Based on Refs. 00, the pair creation operator is defined by 6^ = e^c* 



rn ' 



which generates ( k , f), ( — k , J.) electron pair; and the pair annihilation operator is defined by b m = c m C- m , which 
destroys ( k , t), (— k , |) electron pair. So one can write the Hamiltonian Q as 0: flees = J2m=i £m(n m + n- m )/2 — 



l^Sm ,_, b^bi, where £ m = e m — Ef is the free electron kinetic energy from Fermi surface (ep is the Fermi energy). 
There are two cases for every pair state m: "occupation" and "empty" , which are denoted respectively by spin-up state 

Xi = (^^j , spin-down state xo = ^\ ■ Introducing the so-called spin-analogy corresponding of the pair annihilation 

and the pair creation operators b m and b} n as b m =^ (cr^ — ioy^^j /2 = (7^,6^ => {ox^ +* c 4 m ' > ) /2 = c+, as 

well as the pair number operator n m + ri— m => 1 + ui" 1 ^ , one knows that these fermionic pair operators satisfy the 
commutation algebra: si (2) — {b m , 6^, n m + n_ m — l}, i.e. sl(2) = {er^, <j+, of„}. By use of the above relations, one 
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can express Hbcs i n terms of the spin operators: 

JV N 
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where e m = £ m — V, and the constant term — e m has been ignored Q , which vanishes anyway since we cut off 

m 

symmetrically above and below ef- 

From another way, by use of Bogoliubov transformation and mean field approximation 0, one can obtain the 
diagonal form of BSC-type pairing model Hamiltonian 



1 N 

Hbcs =e s + ^ + A2 ) V2 (7^7 m + l ] - m l-m) (3) 



m— 1 



Here 7^, 7 m are the quasiparticle creation and annihilation operators; e s is the ground state energy of superconducting 
system. It implies that there exists such physical approximation that Hbcs can be diagonalized. Note that pairing 
model Hamiltonian is symmetric and geminate about m and — m, without loss of generality, we are always able to 
assume that it is diagonalized as 



N 
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m=\ a, /3=0 



a/3 (rC^-m + fc) (4a) 
y(»m + n-m) + E E ^ ("™"-™ + ™m™-m) ( 4b ) 



m=l a, 0^0 



where n^j and fp_ m are the number operators of quasiparticle and e a /3 is a symmetric tensor. Now, we use the 
spin-analogy idea to this diagonalized Hamiltonian and then get 

N F' N 

<c a I = ^ + ET L ( 1 + ^ m) )+ E 4,(i + ^ m) )(i + ^ } ) (5) 

m=l ro,(=l,m<i 

where we have used the facts that [(1 + a z )/2] n — (1 + a z )/2 and cr( + 'er'~) = (1 + <r z )/2. Obviously, for the mean 
field approximation, e' rnl f» and E^, = + A 2 ) 1 / 2 . In fact, since the success of mean field approximation, one 
has the reason to believe that e' , almost have no the important effect for the physical spectrum. Of course, if one 
can prove that the diagonalized form of BSC-type pairing model Hamiltonian does not contain the high order terms 
about n m and rt_ m (via. the last terms in eq.©), thus the spectrum in the subspace I of H p is just the exact solution 
of BCS-type pairing model. 

Furthermore, in order to understand the relation between the spectrum of spin-analogy of pairing model and the 
spectrum of pairing model itself, we still need to analyse the structure of the eigenvectors of H p . Note the fact that 
the spin space can be divided into the different subspaces which correspond to the different numbers of spin-up states, 

that is = © S[ © © • ■ • © S^, where the subspace n, i.e Sn, is a subspace with n spin-up states 
|0), we can denote the basis of these subspaces as 

SfW = {|2"», ^ = {11)} (6a) 



o(N) = 1 W \ 
« I iii2---i n / 



>N _Y^2( N - l °A, l a ^l b ,l a = 1,2,- ■■ ,N\ (6b) 




For example, for the subspace I of S^ N \ its basis are ) = (2^ — 2^ ll '),(ix — 1,2,- ■■ ,N). Obviously, 
((2^ — 2( Ar ~ n ))) fc = 5 k 2 N_ 2 < JV - i i) ■ ^ is easy to see that 



rrdiag 

-"bcs 



tf ) )=E' il 6 iljl (7) 
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Then, let us back to our problem, via. to analyse the structure of eigenvectors of H p . From (a^^ai 
{ctx™ 1 + icy 71 ^) — ioy \ > ( m m cq.©, it follows that |1) and [2^) must be H p s two eigenvectors respectively 
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corresponding to the maximum and the minimum eigenvalues. They can be denoted respectively by \vz) and 



Moreover, for the arbitrary basis 



belongs to Sn > H } 



also belongs to S„ because that a + and a 



appear in pairs or do not appear in the various terms of H p . It implies that \s\^ 
,N—1). Therefore 



(JV) 



1,2, 



(JV) 
sub2 



H 



(JV) 
subN 



H, 
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So we can denote the others eigenvectors of H p as v p " J — J^jC") |j^ n ') S with the corresponding 

eigenvalues , and i( n ),j( n ) only take the sequence number of spin basis belonging to the subspace Sn \ for 

example i^(k) = 2 N - 2( N ~ k \ (k = 1,2, • • • , A). In other words, |iW>, |j< n )) € S^. Actually, above statement 
implies that the elements of the diagonalization of H p in the subspace 1 are just E'^ . If one only would like to obtain 
the spectrum of -ffBCS m the physical approximations such as the mean field approximation, it is enough to find the 
eigenvalues of the subspace 1 of H p N ^ and does not need to directly deal with the concrete approximate procedure. 
For our propose, not only we need to express the Hamiltonian Hp' as the direct-sum of the submatrices 0, @, 

For example, when N = 3, the basis of subspace one are 



but also we have to obtain the concrete forms of H 
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|011), 1 101) and 1 1 10} . In the spin space, they are (|011))j = 5 i 2 3_ 2 <3-i) = Sn, (|101)), = <5 ii2 3-2( 3 - 2 ) = <^6j (|H0))j = 

<5j 2 3_2< 3 - 3 ) = Sir, (i = 1, 2, • • • , 8). Thus, denoting H p N ^ is an element in the i-th row and the j-th column H p N \ 
we have 



H 



i?f[4,4] 4 3) [4,6] i?f[4,7] 
£li= | #f [6,4] 4 3) [6,6] ^ 3) [6,7] 
^ 3) [7,4] ^ 3) [7,6] i4 3) [7,7] 
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Likewise, in terms of computer, we can calculate, up to N = 11, H^X has the same structure. That is, all of the 



non-diagonal elements are — V, the diagonal elements are a linear combination of all of e;, (i, 1,2, 
coefficients the i-th being 1/2 and the others are —1/2. That is 



1 N 
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i±3 



-V (i,j = 1,2,- ■■ ,N) 



, N) with the 
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Obviously, it is easy to diagonalize above H^ bl . In our paper we just used this fact to calculate the numerical 
diagonalization of the subspace 1 of H p up to N = 100. Concequentely, the advantage of our method is just to 
transfer an exponentially complicated problem to a polynomial problem. In other words, a task dealing with a 2 100 
dimensional matrix is transferred to a task dealing with 100 dimensional matrix, which can be carried out by the 
usual computer. 

Now, we would like to prove that H^Q^ indeed has the structure in ea. (|10fl in terms of mathematical induction. 



We have obtained that the structure of H^ bl is true when N — 3. Assume that the structure of iJ^Vi ^ s ^id- 
Consider the case N + 1, we can rewrite Hamiltonian as 



(ii) 



where / 2 x2 is a two dimensional identity matrix. 

Denoting pf = 2 N - 2 N ~\{i = 1,2,--- ,N), we have ( 
2(2 N - 2 N ~% it follows that 



,(JV+1) N 



It implies that for a matrix A, 



S ktP N. Obviously, since 2 N+1 - 2 N+1 ~ 



) = K>(°)' |.K 1} ) = |2 (JV+1) -i> 



(A® J 2x2 ) p «+i ]P « + i , (i<iV) 



(12) 



(13) 
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Therefore, the first term Tf 



(JV+l) 



in eq.mj contributes a same N dimensional submatrix as -ff stlbl 



(TV) 



It is easy to see that only the last component of the 2 N -th row and/or the 2 w -th column of Hp N ^ is non zero since 



(m) (0 



<Jx + <J 



(m) m 



(JV+l) 



( Ox 



m) + i4 m) 



g x — i(Ty l> I , (m ^ I) in eq. (J2J) ■ Thus, the contribution of the first term 



.(0 



to the (N + l)-th row and the (N + l)-th column of the submatrix 1 only has a non zero component at the 



1 

last (diagonal element) that is the same as the 2 Ar -th diagonal element of Hp N \ the other components are zero 



Obviously the second term T« N = eAr +1 eri W+1 '/2 in ea. ljlip is a diagonal matrix and for the positions +1 , (i = 
1, 2, • • • , N + 1), except for 2 N+1 — 1 gives a positive ejv+i/2, the others are — £at+i/2. It is easy to see the same 
conclusion in terms of mathematical induction when we consider J 2 x2 <8> T^^" 1 " 1 ^ = T^ N+2 \ 

To join the contributions to N x N submatrix in if^t and the (N + l)-th and/or the (N + l)-th column, and 



(N+l) 



add the contributions from the T- 



(N+l) 



and T, 



(tv+i; 



, we have 



r T {N+i) + t (n+i] 



N+l 



tbl 



[I, I 



m— 1 



(JV+l) 



(JV+l) 



( 

(r 1 (Ar+1) + r 2 (Ar+1) ) [i,jv + i] 



-v, 

r T (N + i )+T ( N+ x)\ [N + 1J] 

V / subl 



= o 



3<JV 



(14a) 

(14b) 
(14c) 



where we have used that the notation denotes the element in the i-th row and the j column of a matrix and 

subscript "i" denotes the i-th component of a vector. 

Now, let us consider the last two terms in ea. lftTj) . We divide three steps to do this. 



Step one: We first prove, in terms of mathematical induction, that t 3 

f N 

E 

\il=l 



i N) \j,2 N ]=4 N \2»,j] = 



»1 
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and 
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(4 W) ) 



j cri"^ satisfies: 
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For 3 spins we have 

/0 1 1 1 0\ 

10 10 10 

1 1 1 

1 1 1 

1 1 1 
1 1 1 
1 1 1 

V 1 1 1 } 

obviously ea. H15|) is correct. Assume that for M spins, ea. (|15f) is also true. Thus, for M + 1 spins, we have 

M+l 



E-i m) 



(16) 
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m— 1 



(M) 



1 ^2x2 



r (M+l) 
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From ( t 



(4 M) ^ 2X2 ) [i,2 



M+l] _ 



and t ri M+1) [z,2 M + 1 l= ( 5 



• 2 M + 1_1, 



it follows that for M + 1 



spins, the first relation in eq.JTEJ is valid. Then, notice that (t ( 3 N) ® 7 2x2 j [pf +1 ,pf +1 ] = 4 N) \pf,pf ] = (i, j < N), 
(4 Af) ®^2x 2 ) [p: 



N+l 9 M+1 



-1]= * 



[2 M+1 — l,Pj +1 ] = 0, the first term in ea-H17|l has no contribution to 

its submatrix 1 . While all of non zero elements of the second term in eq. H17|l is in its two slanting diagonal lines near 
main diagonal line and then has no contribution to its submatrix 1 either. Therefore, the second relation in ea. i|15|) 
for M+l spins is valid. 



Step two: Because of ea.l(T5)l. (j2m=i ^ ' 

Step three: First we prove that 



) [i,j] = 0,(i > j<N). 

/ subl 



/ N 
Vro=l 



/N+l 



f N 

^2 o- lm ■• 

\rn = l 



[2»,i] 
[2 N ,i] 



rN+l 



tN+1 



[M 
[M 




(i 



(N+l) 



,N+1 



<i<2 



N+l 



l,2 iv+i ) 



(18a) 
(18b) 
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Obviously, for 3 spins, we have 




0000000100010100 
0000001000101000 



(19) 



i.e eqs. (I18all8b|) are correct. It must be emphasized that the right side of eas. H18all8b|) corresponds to 2 N+1 — 1 
column in the whole matrix. Suppose that for M spins, eas. <|18all8bl) are correct. It is easy to prove that for M + 1 
spins, this conclusion is correct either since the first relation in ea. l|15[) and the first row of a x being (0, 1). In other 
word, for arbitrary N spins, this conclusion is still valid. Therefore, when we extract the elements of 2 N+1 — 1 row 
and the 2 N+1 - 2 N+1 ~ l columns (i = 1, 2, • • • , N + 1) of H^ N+1) from the third term in ea. (|llf> . we obtain a vector 
whose previous N components are all —V/2, the last component is zero. In the same way, we can prove that when 
the elements of 2 N+1 — 1 column and the 2 N+1 — 2 N+1 ~ l rows (i — 1, 2, • • ■ , N + 1) of Hp N+1 ^ are extracted from the 
third term in ea. l|ll|) . we obtain a vector whose the previous N components are all — V/2 and the last component is 
zero. 

To join the contributions to N x N submatrix in H^J[ and the (N + l)-th and/or the (N + l)-th column from 
T 3 (JV+1) = -V^L, *™*<? +1 >/2, we have 

(r 3 (JV+1) ) [i,j}=0 (i,j<N), , (t 3 (W+1) ) [Y+1,Y+1]=0 (20a) 

V / subl V / subl 

(t^ N+1} ) [i, N + 1] = (t^ N+1) ) [N + 1, i] = -V/2 (t < N) (20b) 

V / subl V / subl 

By use of the same method, it is able to prove that when we extract the elements of 2 N+1 — 1 row and the 
2 n+i_ 2 n+i-i co i umns (j = lj2) . .. ,N+1) or the 2 N+1 - 1 column and the 2 N+1 - 2 N+1 ~ i rows (i = 1, 2, • • • , N+ 1) 

of Hp N+1 ^ from the fourth term in ea. l|ll|) . we obtain a vector whose the previous N components are all — V/2 and 

the last component is zero. That is, the submatrix 1 of T^ N+1 ^ = — VJ2m=i ^v"^ '< J y N+1 ' > /2 is same as one of T^ N+1 \ 
In the last, we add the all of contributions to the submatrix 1 from every term, our proof in terms of mathematical 

induction is finished, i.e the equation (|10l) is correct. 

As an interesting application of our result, we have calculated the energy gap in superconductor [Tol | and aimed at 

the check the precision of quantum simulation in our proposal of quantum simulation of pairing model 0, @ either. 

We believe that there are the other applications in the relevant problems of pairing model. Moreover, it is possible to 

extend our method to the other subspaces as well as the other models of spin systems. 
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